INVERSE BOUNDARY VALUE PROBLEMS FOR THE 
MAGNETIC SCHRODINGER EQUATION 



MIKKO SALO 

Abstract. We survey recent results on inverse boundary value prob- 
lems for the magnetic Schrodinger equation. 



1. Introduction 

The purpose of this article is to give an overview of recent developments 
in inverse boundary value problems for the magnetic Schrodinger equation. 
In these problems, one attempts to find information about the coefficients 
of the equation from measurements made at the boundary. A prototype 
for this class of questions is the Calderon problem [3J, also known as the 
inverse conductivity problem. There one wishes to determine the electrical 
conductivity of a body by measuring electrical currents corresponding to 
voltage potentials at the boundary. Both the theoretical and applied aspects 
of this problem have been extensively studied, see the surveys [M], [35] . 

Inverse problems for Schrodinger equations are closely related to the 
Calderon problem. In fact, most results for the Calderon problem in three 
and higher dimensions are based on a reduction to the Schrodinger equation 

(-A + V> = 0. 

Many other inverse problems, both for scalar equations and systems, have 
been treated using similar reductions. However, often these reductions will 
result in first order perturbations of the Laplacian, which correspond to 
equations of the form 

n f> 
(_A + £A— + V)u = 0. 

3=1 0x * 

The magnetic Schrodinger equation is a convenient model for these equa- 
tions. By studying such first order perturbations, results have been obtained 
in inverse problems for isotropic elasticity [23], [21], [llj . Maxwell equations 
in chiral media [19], Dirac equations [22], [18], Schrodinger equations with 
external Yang-Mills potentials [8], and the Stokes system [12j . 

The magnetic Schrodinger equation is of course important in its own right, 
and it is a natural elliptic equation for which one can study inverse problems. 
It is well known that the inverse boundary value problem for the magnetic 
Schrodinger equation is equivalent to an inverse scattering problem at fixed 
energy, provided that the coefficients are compactly supported. Therefore, 
methods introduced for the inverse boundary value problem are often also 
useful in studying inverse scattering problems. 

l 
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We wish to give a precise formulation of the inverse boundary problem. 
Let f2 C R n be a bounded open set with Lipschitz boundary. We always 
assume n > 3 in this article. The magnetic Schrodinger operator is defined 
as 

H A y = (D + A) 2 + V 

where A = (A±, . . . , A n ) is a vector field (magnetic potential) with compo- 
nents in and V is a function (electric potential) in L°°(Q). The 
coefficients are allowed to be complex valued. Also, D = -V, and V + iA = 
i(D + A) is the magnetic gradient. In nondivergence form, H A y is a first 
order perturbation of the Laplacian, given by 

H Ay = -A + 2A ■ D + V, 

where V = A 2 + D ■ A + V. 

Consider the weak Dirichlet problem 

J Ha,vu =0 in Q, 
\ u = f on dQ. 

By standard methods, Ha.v has a countable set of Dirichlet eigenvalues in fi. 
We will make the standing assumption that is not a Dirichlet eigenvalue. In 
this case the Dirichlet problem has a unique weak solution u = ut G H (CI) 
for any / G H l / 2 {d$). 

The boundary measurements related to H^y are given by the Dirichlet- 
to-Neumann map (DN map) 

Aa,v : / •-»• (V + iA)uf ■ u\ 9n . 

Here v is the outer unit normal of d£l. Thus, for any boundary value /, one 
can measure the magnetic normal derivative of Uf at the boundary. More 
precisely, we define A^y using the weak formulation 

(A Ay f\g) dn = ((V + iA)u f \(V + iA)e g ) + (Vu f \e g ), f , g e H 1 / 2 (dQ) , (1) 
where e g is any function in H 1 ^) with e g \dQ = g, and where we write 

(u\v)= / u-vdx, (f\g)asi= / fgdS. 
Jn Jan 

Here dS is the Euclidean surface measure on <9f2. It follows from Q that 

A A ,v is a bounded map H l / 2 (dQ) -> fl" -1 / 2 ^). 

An important property of the magnetic Schrodinger operator is gauge 

invariance. This means that in the weak sense 

e-^H Ay (e^u) = H A+S7tp>v u (2) 

if (p e W 1,0O (n) and u G i? x (0). If additionally <p\an = 0, the boundary 
measurements are preserved: 

Ayl+V V ,V = A A y. (3) 

This gauge transformation replaces the magnetic potential A by a gauge 
equivalent potential A + Vtp. However, the magnetic field dA is preserved, 
where dA = d(Y^=i Aj(x) dxj) is the exterior derivative of A interpreted as 
a 1-form. If n = 3, then dA corresponds to V x A. 

We may now formulate the basic inverse problem considered in this article. 
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Inverse problem for the magnetic Schrodinger equa- 
tion: Given the boundary measurements Aa,v, is it possible 
to determine the magnetic field dA and the electric potential 
V in n?. 

As with other inverse problems, there are different aspects which can be 
studied. We will concentrate on the following topics. 

(1) Boundary determination: given A^y, determine the values of A 
and V on d£l. 

(2) Uniqueness: if h-A\,Vi = ^-A 2 ,V2-> show that dA\ = dAi and V\ = V%. 

(3) Reconstruction: given Aa y, reconstruct the values of dA and V . 

(4) Stability: if A_Aj_ vj and Aa 2 ,v 2 are close, show that dAi and dA2 
(and V\ and V2) also have to be close. 

(5) Partial data: if A^yjr = A J 4 2i v - 2 |r for a subset V C c?S7, show that 
dA\ = dA 2 and V\ = V 2 . 

The inverse problem for the magnetic Schrodinger equation has been stud- 
ied by several authors. One of the first results is due to Sun [29], who showed 
uniqueness for magnetic potentials in W 2,co satisfying a smallness condition. 
The proof was based on complex geometrical optics (CGO) solutions, which 
were introduced by Sylvester and Uhlmann [30] in the case ^4 = 0. If A is 
nonzero and large, then the construction of CGO solutions is more difficult, 
and was achieved by Nakamura and Uhlmann [23], [24] using a pseudo- 
differential conjugation method. With this method, Nakamura, Sun, and 
Uhlmann [21] showed uniqueness for smooth coefficients and domains with- 
out any smallness conditions, and they also gave a boundary determination 
result. The uniqueness result was extended to C 1 magnetic potentials by 
Tolmasky [32], and to some less regular but small potentials by Panchenko 

The fixed energy inverse scattering problem for magnetic Schrodinger 
operators has been studied by Henkin and Novikov [25] for small A, and 
by Eskin and Ralston [S] without smallness assumptions. There is a close 
relation between the inverse boundary value problem and the fixed energy 
inverse scattering problem, and in fact for compactly supported potentials 
the two problems are equivalent. See Eskin and Ralston [TO] for a proof of 
this in the magnetic case. We will not consider scattering problems in this 
article. 

Recently, there have been many developments in the above questions. 
Brown and Salo [2] proved a boundary determination result for continu- 
ous magnetic potentials in C 1 domains. Interior uniqueness was shown to 
hold for Dini continuous magnetic potentials by Salo [27] . Salo [28] gave 
a reconstruction algorithm, and stability results were proved in Tzou [33]. 
A partial data result was given by Dos Santos Ferreira, Kenig, Sjostrand, 
and Uhlmann [6] for C 2 magnetic potentials, following the result for A = 
by Kenig, Sjostrand, and Uhlmann [15] . These works also introduced a 
new method for constructing CGO solutions, based on Carleman estimates 
and a convexification idea. The partial data result was extended to Holder 
continuous potentials by Knudsen and Salo [16]. 
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We will describe these recent developments below. The article is orga- 
nized as follows. Section [2] explains the boundary determination result in [2\. 
In Section [3l we discuss the construction of CGO solutions using pseudodif- 
ferential conjugation, following [28]. Section H] is devoted to the uniqueness, 
reconstruction, and stability results given in [27], [28], [33]. In Section [5] we 
describe the construction of CGO solutions based on Carleman estimates 
following |6j and [16], and the final Section considers the partial data 
results given in these papers. 

Notation. All functions will be complex valued, unless stated otherwise. 
We let C k (Q) be the set of functions whose partial derivatives up to order k 
are continuous in O, and C k+£ (Q) consists of the functions in C k (Q) whose 
partial derivatives of order k are e-Holder continuous in Q. The set of 
compactly supported continuous functions in R n is denoted by C c (R n ), 
similarly C£°(R n ) and L?(R n ). The L 2 -based and L p -based Sobolev spaces 
are denoted by H s and W k,p , respectively, and Hq(Q) is the set of functions 
in with zero boundary values. Functions with values in a Banach 

space X are denoted by C k (Q;X) etc. 

2. Boundary determination 

In the context of the Calderon problem, if the domain and conductivity 
are smooth, it was shown in [3TJ that the DN map is a pseudodifferential 
operator, and the Taylor series of the conductivity at the boundary can be 
read off from the symbol. The same philosophy applies to many inverse 
problems. For the magnetic Schrodinger operator, this method was used 
in [21] to show that if everything is smooth, then the Taylor series of the 
tangential components of A and of V at the boundary can be determined 
from Aav (there is a small mistake in [21\ Theorem D], but this is not hard 
to fix). The tangential components of A on d£l are 

^4tan = A — (A - v)v. (4) 

The gauge invariance ([3]) shows that one can only expect to recover A ta n 
from Aa,v- 

If the domain and coefficients are not smooth, then pseudodifferential 
theory is not available. However, the symbol of a pseudodifferential operator 
can be obtained by testing against oscillatory functions, and this idea works 
also in the nonsmooth case. Early results are given in [T7]. For the Calderon 
problem in Lipschitz domains the method was used in [TJ , and the magnetic 
Schrodinger equation in C 2 domains was considered in [27]. The following 
theorem is proved in [2]. 

Theorem 2.1. Let O be a C 1 domain and n > 3. Let A £ C(U; C n ) and 
V £ L°°(J7). Then Aa,v determines A tw on dil. 

More precisely, let xq G dfl and let a be a unit tangent vector to dQ at 
xo- In [2] one constructs a family of oscillating functions Jm G L 2 (dQ), with 
/m supported in B(xq, 1/M) n dQ, such that 

lim ({Aa,v ~ A ,o)/m|/m)^ = a • A(x ). (5) 

M—>oo 
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The method is local: one only needs to know A^yf near xq for functions / 
supported near xq. The method also gives the stability estimate 

11(^1 ~ A 2 )ta.n\\L°°(dn) < V2\\A Al y 1 - Aa 2 ,V 2 ||l 2 (9Q)^L 2 (9Q) > 

if Aj are continuous in £1 and Vj are L°°. 

We now describe some ideas in the proof. Let £1 be given by a defining 
function p G C 1 (R n ;R), so that = {p > 0} and 9^ = {p = 0}, and 
V p 7^ on dfl. Let xo and a be as above, and normalize p so that Vp(xo) = 
— ^ (xo). We also choose coordinates so that xq = and z^(xo) = — e n , the 
nth coordinate vector. 

Suppose uj is a modulus of continuity for Vp, and let 77 G C£°(R n ) be a 
function supported in B(0, 1/2) and satisfying J Rn -i rj(x' ,0) 2 dx' = 1. Put 
Vm(x) = r](M(x', p(x))) where M is a large parameter, and define 

u (x) = r jM (x)e N ^ a - x -P^ , x G R", 

n — 1 

/m (x) = M— u (x), x G dtt, 

where N is a parameter larger than M, defined by 

M" 1 o;(M- 1 ) = N" 1 . 

The idea is that uq is an approximate solution to H^yu = in $7, concen- 
trated near the boundary point and oscillating at the boundary. 
From (pQ) we obtain 

((A A y-A . )f M \f M ) d n = M n - 1 [ {iA-(uVv-vVu) + {A 2 + V)uv)dx 

Jn 

(6) 

where u G solves Hyyu = with u = uq on <9f2, and v G 

solves —Aw = with ?j = uq on <9f2. We write u = uq + ui and w = uo + v 1 
where Ui,«i G ffg(O). Since no was chosen as an approximate solution of 
H^yu = (and — An = 0), u\ and v\ will be error terms smaller than uq, 
which solve 

Ha,vUi = -H A yu in Q, (7) 
-Avi = Au in tt. (8) 

Using the explicit form of uq, the main term in © will satisfy 

lim M™ -1 / M • (n Vn - uoVuo dx = ,4(0) • a. 

Thus, ([5]) and Theorem 12.11 will follow if one can show that all the terms in 
([6]) involving u\ or v\ will be o(l) as M —* 00. 

Since ui,t>i are i?o(£2) solutions of (|7|), ([8]) where the right hand sides 
are explicit, we get good bounds for ||Vui||^2(q) and ||Vt> lH^am) from weak 
solution estimates. The problem is to estimate ||iti||i2m) and Huill/^^)- We 
will use Hardy's inequality, 

IIW^IIl 2 ^) < C\\Vw\\ L 2 {n) , w G Hq(SI), 

where <5(x) = dist(x,50), and also the estimate 

\\SVw\\ L 2( n ) < C\\w\\ L 2(q), w G L 2 {Vt) harmonic, 
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which follows from L 2 boundedness of the Hardy-Littlewood maximal func- 
tion (see [27j Lemma 5.5]). Using these two estimates in ([6]), it will be 
enough to obtain a bound for ||ui||£2(n). 

The main technical point of [2] is to prove for C 1 domains the estimate 

1— n 1 , , 

lh|| i2(n) < CM— N~ 5. (9) 

This will allow to finish the proof of Theorem l2.ll To prove ([9]) we first note 
that the Poincare inequality and the gradient bound for v\ yield a good 
estimate for \\v\ ||L 2 (nns(o,iooM- 1 ))- In fi \ -6(0, 100M -1 ) we know that v\ 
is a harmonic function vanishing on d£l. Now we use the fact that there 
is another such function for which good bounds are known, namely the 
Green function G(x, y) for the Laplacian in f2, with singularity at a suitable 
point in B(0, lOOAf -1 ). The pointwise values of v\ can be compared with 
G, using the boundary Harnack principle of Jerison and Kenig [13J. In 
high dimensions, one gets a bound for ||^i Hi 2 (O\S(0,i00JW— !)) from the easy 
estimate 

\G(x,y)\ <C\x-y\ 2 ~ n . 

In lower dimensions one needs to use additional vanishing of G at the bound- 
ary. We refer to [2] for the details. 

3. CGO SOLUTIONS AND PSEUDODIFFERENTIAL CONJUGATION 

In [30], CGO solutions were introduced for the equation Hoyu = 0. They 
are of the form 

u = e lp - x (l + r p ) (10) 

where p £ C n is a complex vector satisfying p ■ p = 0, and r p is a remainder 
term which is small when \p\ is large. Thus, for \p\ large the solution u 
resembles the harmonic exponential e lp ' x , and one can use these exponentials 
to recover the potential V from boundary measurements Aoy. The method 
of CGO solutions is central in most results on inverse boundary problems 
for elliptic equations. 

We wish to discuss CGO solutions to the magnetic Schrodinger equation 
Ha,vu = 0, or equivalently to (—A + 2A ■ D + V)u = 0. This is more 
difficult than for Hoy. The first construction was achieved by Nakamura 
and Uhlmann [23], [21], who introduced a pseudodifferential conjugation 
argument which essentially removes A. Stronger forms of this construction 
are now available. In |32j the argument was extended to C 2 / 3+e magnetic 
potentials, and in [27] to just continuous A. These results only yield solutions 
in a bounded domain, and there is no uniqueness for the solutions. In [28] 
the construction was extended to yield global unique CGO solutions. 

We will follow [28] . First we introduce weighted 1? spaces L|(R n ) with 
norm \\{x) s f\\i2^n^, where (x) = (1 + l^ 2 ) 1 / 2 , and weighted Sobolev spaces 
#£(R n ) with norm — Xyla^fcll^/ll-Ll' The following is the main 

result in this section. It shows existence and uniqueness of CGO solutions 
for large \p\ with known asymptotics as \p\ — > oo. 



INVERSE PROBLEMS FOR THE MAGNETIC SCHRODINGER EQUATION 



7 



Theorem 3.1. Suppose A G C c (R n ; C n ) and V G L£°(R n ). Let -1< (5 < 
0. If p G C n satisfies p ■ p = and if |/o| is large enough, then the equation 

(-A + 2A ■ D + V)u = in R n 

has a unique CGO solution of the form (|10p with r p G H*. Further, if p has 
the form p = j:{yi +1V2) where u\, v% G R n are orthogonal unit vectors, and 

If 1 

< t>{x) = -— / -^Oi + ^2) -^(s - yi ^1 -2/2^2) dyidy 2 , (11) 

then one has the asymptotics 1 + r p = a p + f p where 

a p — > e i?i pointwise, ||ap||i°o = 0(1), /i||Va p ||L°o = o(l), (12) 
||fp|| L 2 +/i||Vf p || L 2 = o(l), (13) 

as h — > 0. 

To discuss the proof of Theorem[37TJ we introduce the conjugated operator 

e- ip - x H Ay e ip - x = P P + 2A-D p + V 

where P p = -A + 2p ■ D, D p = D + p, and V = A 2 + D- A + V. Inserting 
(flO|) in Ha,vu = 0, we see that the construction of CGO solutions reduces 
to solving 

(P p + 2A-D p + V)r p = f (14) 

where / = -2A ■ p - V. If A = V = the solvability of dHJ) follows from 
the fundamental estimates of Sylvester and Uhlmann. 

Proposition 3.2. [30] Let -1 < 5 < 0, and let p G C n with p - p = X and 
|p| > 1. Then for any / G Lg +1 the equation P p u = / has a unique solution 
u G L^. The solution operator G p : / 1— » n satisfies 

1^111^/11^ + ||VG p /|| L2 <C||/|| L2+i . 

If ||^4||x,oa is small then we may solve (I14p by trying r p = Gp^p and using 
Neumann series. If A is large this will not work, since then I + 2A ■ D p G p 
is not a small perturbation of identity on L|m. 

If A is large, one might try to use a gauge transformation as in ([2j) to 
remove A However, if A is not a gradient, this will not remove all of A. 
The idea of Nakamura and Uhlmann was to replace the exponentials e lip in 
([2]) by more general pseudo differential operators. This may be thought of 
as a pseudodifferential gauge transformation, which essentially removes A. 
The rest of the proof then proceeds by a perturbation argument as in the 
case where A is small, using the norm estimates for G p and Neumann series. 

It will be convenient to switch to semiclassical notation. This amounts 
to writing h = r^r and to considering a small parameter h instead of a large 
parameter \p\. However, we may then use the well-established machinery of 
semiclassical pseudodifferential calculus, which automatically keeps track of 
the dependence on h in norm estimates, and simplifies some proofs since the 
parameter h is scaled away when passing to symbols. 

We will use the usual semiclassical symbol classes, see [5]. 
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Definition. If < a < 1/2 and m G R, we let S™ be the space of all 
functions c(x,£) = c(x,£;/i) where x, £ G R n and h G (0,ho], ho < 1, such 
that c is smooth in x and £ and 

\d^c(x,0\<C aP h-^\(O m 
for all a, f3. If c G S 1 ™ we define an operator C = Op h (c) by 

C/(x) = (27T)-" jj^/ {X ~ y)<c {^i h t) f(v)^dy. 

Semiclassical quantization means that the symbol £j corresponds to the 
operator hDj. Note that we have used Weyl quantization. The operators 
have a calculus: if c £ S® then C is bounded on weighted spaces L 2 with 
norm independent of h, and if c G S 1 ™ and d G 5™' then CD is an operator 
with symbol in 5™+ m . 

To manage the nonsmooth coefficients, we introduce the standard mollifier 
Xs(x) = 5~ n x(x/5) where x G C^°(R n ), < x < 1; an d J xdx = 1. Fix do 
with < do < 1/2, and consider the /i-dependent smooth approximations 

A« = A*x* (15) 

with the specific choice (5 = . 

We may write the equation (|14p in semiclassical notation as 



(P + hQ + h 2 V)r p = h 2 f (16) 

where P = -h 2 A + 2p-hD and Q = 2A • (/iP + p), and p = hp so |/5| = 1. 
We split Q as Q = Q tt + Q b where Q tt = 2A» • (/iP + p). The following result 
states that the smooth first order term Q$ can be conjugated into a term 
which has order less than 1. 

Proposition 3.3. Given a with uq < a < 1/2, there exist c, c, s G S 1 ^ and 
(5 > such that 

(P + hQ*)C = CP + h 1+s (x)- 1 S. (17) 

Further, C and C are elliptic, in the sense that c and c are nonvanishing for 
small h. 

To prove Proposition 13.31 one takes some c G 5° and commutes P + /iQ" 
and C using the symbol calculus. The result is 

(P + hQ$)C = CP + hO-p h (-H p c + q 9 c) + lower order terms. 

Here p(f ) = £ 2 + 2/5 • £ is the symbol of P, g fl = 2A» • (f + iVp) - hD ■ A$, 
and Hp = V^p • V x — V x p ■ = 2(£ + p) • V x is the Hamilton vector field. 
We want to choose c so that the first order term disappears, i.e. 

-PpC + ^c^O. (18) 
i 

Trying c = e lip , this equation becomes 

H p ip « -<?». 

Near p _1 (0) this looks like a 9-equation in some variables, but degenerates 
away from p _1 (0). However, the operator P is elliptic away from p _1 (0), 
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and it will be sufficient to find <p near p _1 (0). Inserting a suitable cutoff and 
applying the Cauchy transform, we obtain a function ip satisfying 

|fi£9fv(*,0l < c^h-^^ix)^- 1 . 

The problem is that the ^-derivatives of ip grow in x, so ip is not in S® . This 
can be fixed by introducing another cutoff \ £ C^°(R n ) with X = 1 near 0, 
and by taking 

c = jx{h 6 x)tp 

where 6 > is chosen so that a := (Tq+9 < 1/2. With this choice c G 5°, and 
we have traded off the undesired x-growth into slightly worse /i-dependence 
in the symbol. This method was used in the context of energy estimates for 
nonlinear Schrddinger equations by Kenig, Ponce, and Vega [14] . 

We have arrived at the correct conjugating operator C, and it is possible 
to finish the proof of Proposition 13.31 Theorem 13.11 follows by combining the 
pseudodifferential conjugation with a perturbation argument where small 
errors are inverted by Neumann series. The asymptotics for CGO solutions 
may be computed either from the asymptotics of C, or by guessing the right 
asymptotics and using the uniqueness part of Theorem 13.11 The details are 
given in [2"g] . 

4. Uniqueness, reconstruction, and stability 

The CGO solutions constructed in Theorem 13.11 can be used to solve the 
uniqueness, reconstruction, and stability questions for the inverse problem 
for H^y. We begin with reconstruction. The next result was proved in [28j . 
We state it here in a simplified form to avoid technicalities. 

Theorem 4.1. Let Q, C R n , n > 3, be a bounded simply connected domain 
with smooth boundary. Suppose A G C c 1+£ (0;C n ) and V G £°°(fi;C). 
Then one may reconstruct dA and V from Kav 

In the case ^4 = this result is due to Nachman [20], who used CGO 
solutions and ideas from scattering theory to reconstruct V from A-^y. The 
fact that one has unique global CGO solutions for Hoy is crucial for the 
algorithm. Theorem 13.11 provides such solutions for H^y, and using these 
we can extend Nachman's algorithm to the magnetic case. 

If p G C n satisfies p ■ p = and \p\ is large, we denote by u p the unique 
CGO solution to Ha,vu = given by Theorem 13. 11 The first step is to show 
that Up may be characterized in terms of an integral equation on d£l. More 
precisely, u p \qq is the unique solution in H 3 / 2 (d£l) to the equation 

(I + S p (A A y-A 0fi ))f = e i P x ondfl, (19) 

where S p is the single layer potential 

S p f(x)= [ G p (x,y)f(y)dS(y), f G H^dSt), x G dQ, 
Jan 

and G p is the p-dependent Green function 
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Note that if p = this would be the usual Green function for —A. 

The point is that if we know the boundary measurements A.a,v, and if 
p is fixed, then we know the quantities in (fT9|) and can determine u p \qq, 
as the unique solution of this equation. Then we also know the scattering 
transform tA,v, defined by 

tA,v(Z,p) = ((Aa,v - \-\tf)(u P \an),e ix <^) dn , 

where £ G R n is a vector, and p G C n satisfies p ■ p = 0, Re p _L £, Im p _L £, 
and |p | is large. 

The expression i^y is chosen so that in the limit \p\ —> oo, we obtain cer- 
tain nonlinear Fourier transforms of components of dA. This motivates the 
name "scattering transform", since similar phenomena occur in scattering 
theory. In fact, using the definition (fT]) of Aa,V and the asymptotics for u p 
given in Theorem 13.11 one sees that 

lim/iUy(£,r(>i + " / 2)) = 2 / e- ix ^e i ^(u 1 +iv 2 )- Adx (20) 

h-^Q ' h J 

where £, v\, v 2 are orthogonal and \v\\ = \v 2 \ = 1, and (/> is as in (jlip . If one 
replaces e J<?i by 1 this would just be the Fourier transform of a component 
of dA, so this indeed looks like a nonlinear Fourier transform of dA. 

However, due to a result in [S], this purported nonlinear Fourier transform 
is just the usual one: integration by parts shows 

J e^VVi + -Adx = J e~ ix ^{vi + iu 2 ) ■ Adx. 

Thus, given A.a,v, we have recovered the last expression for different or- 
thogonal triplets £,z^i,^2 where \v\\ = | z^2 1 = 1- These are just the Fourier 
transforms of components of dA, and we have reconstructed the magnetic 
field. 

After finding dA we want to reconstruct V. First we construct some 

C l+e 

magnetic potential A with dA = dA and A\g^ = 0. Since f2 is sim- 
ply connected, we know that A = A + where p\qq = 0, possibly after 
substracting a constant. By gauge equivalence, we know A^o = A^ Q . For 
£ G R n , define another scattering transform 

tA,v(0 = ((^a,v - AA,o){up\dn),Vp\ d n) d a 

where p, p are complex vectors chosen in a certain way so that p+p = and 
\p\ = \p\ = V^> and Up and vp are suitable CGO solutions whose boundary 
values are obtained as solutions of integral equations. It follows that 

lim t A ,v(0 = [ e- ix< V(x) dx, 

h—*ca ' J 

and the electric potential is determined. 

If one is only interested in uniqueness, the regularity conditions can be 
relaxed somewhat. We say that / is Dini continuous, written / G C Dim (r2), 
if \f(x) — f(y)\ < uj(\x — y\) for some modulus of continuity u, which is a 
continuous nondecreasing function satisfying 

uj(0) = 0, f ^-dt<oo. 
Jo t 
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The following result is from |27j . 

Theorem 4.2. Let £1 C R n be a bounded open set with C 1 boundary, 
and n > 3. Suppose A 1 ,A 2 G C Dini (Tl;C n ) and VI, V 2 G L°°(^;C). Let 
Aai,Vi = Ay^y^. Then dA\ = dA 2 and Vi = V%. 

To prove this, one first uses boundary determination (Theorem 12, ip and 
gauge invariance to get A% = A 2 on dQ, and then the vector fields are 
extended to a larger ball B so that they coincide outside 0, . The DN maps 
on dB coincide since the original DN maps do, and after another gauge 
transformation we may assume that Aj G C(B;C n ) and D ■ Aj = (this 
is where Dini continuity is required). The construction of CGO solutions 
in Theorem 13.11 can now be applied, and we may use the CGO solutions 
and arguments similar to those given in the reconstruction to show that 
dAi = dA 2 and V\ = V 2 . 

We next consider stability results, which state that two sets of coefficients 
should be close if the DN maps are. It is known that many inverse problems 
are ill-posed, and for instance the Calderon problem only enjoys a logarith- 
mic stability estimate. Logarithmic stability also holds for the magnetic 
Schrodinger equation. This was proved in [33J along with a log log stability 
estimate for partial data. We state the full data result in a slightly simplified 
form. 

Theorem 4.3. Let f2 C R n be a bounded open set with smooth boundary, 
and n > 3. Suppose Aj G C'°°(0;R n ) and Vj G L°°(0), and suppose 
ll^H^a.oo < M, \\Vj\\ L oo < M for j = 1,2. There exist C > and e > 
such that 

IKAx - A 2 )|| H _i + \\Vt - V 2 \\ H -x < C(\\A Al , Vl - Aa 2 ,v 2 || 1/2 

+ |log ||A Al! Vi - A A2 ,v 2 \\\~ £ ) 

where || • || = || • \\ H ^{d€i)-^H-y^{dQ) on the ri g ht - 

A main step in the proof is to ensure that whenever the coefficients of 
Ha,v satisfy a priori bounds, such as 

P|| W 2,oo < M, [|V[| L oo(n) < M , A and V supported in 5(0, M) 

where M > 0, then the construction of CGO solutions in Theorem 13 . 1 1 goes 
through with constants only depending on M. This means that there are 
Cm > 0, em > 0, so that the solutions exist for \p\ > Cm and the remainder 
f p satisfies 

\p\¥ P \\q + ¥ P \\Ll <C M \p\ l ~ £M - 

Given this, one may insert suitable CGO solutions in the weak definition of 
Aai,Vi — A J 4 2j y 2 , and the Fourier transform of components of d(A\ — A 2 ) will 
be controlled by quantities depending on HA^y-j — A^ 2 y 2 || and M. For the 
electric potentials, one needs additional gauge transformations to control 
\\A\ — A 2 \\l^ by \\d(Ai — A 2 )\\lp where p > n. 
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5. CGO SOLUTIONS AND CARLEMAN ESTIMATES 

Recently, a new method was introduced for constructing CGO solutions 
to the magnetic Schrodinger equation. The method is based on Carleman 
estimates, and it avoids the use of pseudo differential operators. We will 
describe the main ideas of the construction here. In the next section, we 
discuss how these ideas can be used to prove partial data results for Ha,v- 

If one is working in a bounded domain f2 C R n , the Sylvester-Uhlmann 
estimates in Proposition 13.21 may be written as 

h\\u\\ + h\\hDu\\ < C\\eV /h {-h 2 k)e-v /h u\\, u £ C™(Q), (21) 

where || • || = || • ^Pi. x ) = ot ■ x with a G R n a unit vector, and h > 

is small. Such weighted norm estimates with small parameter are called 
Carleman estimates. Now (|21|) is sufficient for constructing CGO solutions 
of the form (|l(jp to Hoyu = in $7. This point of view was adopted by 
Bukhgeim and Uhlmann [3] , who used (|2*T|) to give a partial data result for 
the operator Hoy. 

Kenig, Sjostrand, and Uhlmann [15] introduced new techniques in inverse 
boundary problems using the Carleman estimate approach. They consider 
more general limiting Carleman weights ip for which estimates like (I2ip can 
be proved, both for ip and — ip. The main examples are the linear weight 
<p(x) = a ■ x, and the logarithmic weight 

<p(x) = log \x — xo|, xo ^ fl (22) 

Approximate CGO solutions to H^yu = are obtained from a WKB con- 
struction for the conjugated operator e^^h^Hoye - ^^ , and the approxi- 
mate solutions are converted into exact CGO solutions by solving an in- 
homogeneous equation. The last step is made possible by the Carleman 
estimate, which also yields decay of the error term in the I? norm. 

Proving the estimate (|21j) usually requires a convexity condition on the 
weight if. Since the estimate also needs to hold for — cp, the limiting Carle- 
man weights <p will only satisfy a degenerate convexity condition. In [15J a 
convexification argument was used to prove pip : one replaces the limiting 
Carleman weight ip by a convexified weight ip s , for instance 

2 

<Pe = <P + h^. (23) 

One then proves the Carleman estimate for ip e , and uses this to deduce the 
estimate for the original weight <p. 

The preceding discussion applied to CGO solutions of H^yu = 0. For 
the magnetic operator Ha,v where A is large, (|2ip does not immediately 
provide CGO solutions, since the first order term is not a small perturbation. 
However, if one uses the convexified weight ip e then (|21|) can be proved also in 
the magnetic case, and this will imply the Carleman estimate for H^y with 
the original weight cp. In other words, convexification makes the procedure 
more robust, so that one can handle first order perturbations. This method 
was used in [6] to construct CGO solutions to Ha,vu = 0, and to give a 
partial data result for the magnetic Schrodinger equation. 
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The results in [6] assume that the magnetic potential is C 2 . In [16], the re- 
sults were extended to Holder continuous potentials. The problem in passing 
from C 2 to C e magnetic potentials is that the regularity of CGO solutions 
deteriorates. In particular, the solutions will not be in H 1 ^). However, 
by combining Carleman estimates with the pseudodifferential conjugation 
argument, we obtain CGO solutions in H 1 ^) even for Holder coefficients. 
One of the main points of [16J is an extension of the pseudodifferential con- 
jugation argument to logarithmic Carleman weights. 

In the rest of this section, we give some more details on the construction 
of CGO solutions using Carleman estimates. The main result will be the 
following. This was proved for A G C 2 in [6], and for A G C £ in [16j . 

Theorem 5.1. Let A G C £ (Tl; C n ) and V G L n (0). Let ip and ip be defined 
by (|22|) and (|31|) . respectively. Then for h small there is an H 1 ^) solution 
u = e~h^ +l ^ (a + r) of the equation H^yu = in 17, where a is given by 
(|32|) . One has the norm estimates 

\\d a a\\ Laa{n) = 0(h-°\«\), 

+ \\hVr\\ L 2 {n) = 0(h° £ ) 

where a > is small. 

We follow [E] and [6] and assume first that A G C 2 . Let tt,n C R n be 
bounded domains with 0, C f2, and let ip G C°°(r2;R) with V<£ 7^ in Q. 
We want to find conditions on ip so that (|2ip holds for ip and —p. Consider 
the conjugated operator 

P^ = e v/h (-h 2 A)e-^ h . 

We would like to show H-P^itH > ch\\u\\ when u G C^°(Q). To do this, we 
write = A + iB where A, B are self-adjoint, and note that 

||iVH 2 = {P v u\P v u) = \\Au\\ 2 + \\Bu\\ 2 + (i[A,B)u\u), 

where [A, B] = AB — BA is the commutator. To get a bound from below, 
one would like a positivity condition for i[A, B\. The principal symbol (in 
semiclassical Weyl quantization, see Section [3]) of i[A, B] is h{a, b}, where 
{a, b} = V^a • V ' x b — V x a ■ Vgfr is the Poisson bracket, and 

a(x,0 =f-(Vp) 2 , b{x,i) = 2Vp-i (xe^GR"). 

The positivity condition for _B] at symbol level is 

{a, 6} > when a = b = 0. 

Since [A, B] changes sign when p is replaced by —p, one also needs {a, b} < 
when a = b = 0. This leads to the following condition. 

Definition. Let p be as above. We say that ip is a limiting Carleman weight 
(for the Laplacian) if 

{a, b} = when a = b = 0. (24) 

The condition (|24p is used in many ways in the construction of CGO 
solutions and determination of coefficients. However, the condition is also 
quite restrictive. 
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Given (|24p . we want to prove a Carleman estimate. We first convexify 
the weight and replace ip by <p £ , defined as in (f23|) . If A £ and B £ are the 
operators corresponding to A and B, one still has 

\\P Vs uf = \\ A euf + \\B £ uf + (i[A £ , B £ ]u\u). (25) 
Now, if e is small enough, one can prove the lower bound 

{i[A £ ,B £ ]u\u) > ^-|M| 2 - ^\\A £ u\\ 2 - ^\\B £ u\\ 2 . 
Note that the first term is positive. This yields the estimate 
4=(IMI + \\hDu\\) < C\\e^ /h (-h 2 A)e-^ /h u\\. 

By choosing e small enough, we can replace — /i 2 A by h 2 H A y. Finally, since 
e tpc/h _ e <p/h m w h ere m i s bounded from above and below in we arrive 
at the Carleman estimate 



u 



+ \\hDu\\) < C\e 9th tfH Ay e-vl h u\, u G C c °°(0). (26) 



Shifting the estimate to a different Sobolev index and using the Hahn- 
Banach theorem, one may show that this a priori estimate implies solvability 
for an inhomogeneous equation. 

Proposition 5.2. [6] Let ip be a limiting Carleman weight in f2, let A G 
C 1 ^), and let V G L°°(0). If ft, is small enough, then for any / G L 2 (Q), 
the equation 

e^ h h 2 HA,ve~ v/h r = f in O 
has a solution r G i^ 1 (r2), which satisfies h\\r\\ + /i||/iDr|| < C||/||. 

Proposition 15.21 will be used to convert an approximate CGO solution 
to an exact one. Given a limiting Carleman weight (p, the CGO solution 
u = e~ Lp l h v will satisfy 

e v/h h 2 H A ye-* /h v = 0. (27) 
We find v by a WKB construction in the form 

v = e ~s^(a + r) 

where ip G C°°(f2;R.) is a phase function, a is an amplitude, and r is a 
correction term. Inserting the WKB ansatz in (I27h , and grouping the terms 
in powers of h, one obtains 

(-(Vp) 2 + ih[Vp o D + D o Vp + 2Vp • A] + h 2 H Ay )(a + r) = 0, 

where p = 92 + This will be satisfied if ip, a, and r are chosen so that 

(Vp) 2 = 0, (28) 
(Vpo D + D oVp + 2Vp- A)a = 0, (29) 

e i el h h 2 H Ay e-^ h (e- i t r ) = —h 2 e~ i ^hH A ya. (30) 

The first equation may be written as 

(Vip) 2 = {V<p) 2 , V<p-Vip = 0. 
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This is an eikonal equation for tp. If tp is known, (|29p is a transport equation 
for a, and when tp and a have been found, the last equation for r can be 
solved by Proposition 15.21 

We now choose <p as in ([22]) , and indicate how to solve ([28]) - ([29]) following 
[15] . [6]. It is known that distance functions are solutions of the eikonal 
equation, and one may take 

ip(x) = distcn-i ( ^-,uj), (31) 

V \x — xq\ J 

where uj G S n ~ l is chosen so that tp is smooth near £2. If one chooses 
coordinates so that xq = 0, £1 lies in {x n > 0}, u = e%, and x = {x\,r9) 
with r > and 6 G S 1 ™ -1 , then 

(p + iip = log z 

where z = xi + ir is a complex variable. In these new coordinates, the 
transport equation ([29]) becomes 



% Tl — 2 

(d- z + ~(e 1 + ie r ) ■ A(z,8) - ^ _ )a(z, 6) = 

where e r = (0, 9). This has a solution 

a(z,6) = (z - z)^e i$ (32) 

where $(z, 0) = — ^37 1 ((ei+ie r )-A(z, 0)) is obtained as a Cauchy transform. 
This ends the proof of Theorem 15.11 in the case A £ C 2 . 

If t4 is only C e , the main step is to prove Proposition 15.21 under this 
assumption. We follow [16] and write, as in (]T6l) . 

e^ h h 2 H Ay e-^ h = P + hQ + h 2 V, 

where P = e^ h (-h 2 A)e-^/ h and Q = e^ h (2A ■ hD)e~^ h . The free equa- 
tion Pr = f can be solved by Proposition 15.21 and one wishes to use 
a perturbation argument to solve the full equation. Again, the first or- 
der term Q is not a small perturbation. We split Q = Q" + Q where 
q$ = e ^/ h [2A^ ■ hD)e~ ip l h , with A^ a /i-dependent smooth approximation as 
in ()15p . with 5 = h°~ and a > small. Then will be a small perturbation 
which can be inverted by Neumann series, and the smooth part can be 
conjugated away using pseudo differential operators as in Proposition 13.31 

Proposition 5.3. There exist c, c, r £ S® so that 

(P + hQ^C = CP + h 2 - 2a R in ft. 

Both C and C are elliptic, in the sense that c and c are nonvanishing for 
small h. 

The proof is similar to that of Proposition 13.31 Here one does not need 
global solutions, so it is enough to use cutoffs to restrict to bounded domains. 
However, the transport equation (TT8]) does not have constant coefficients 
anymore, and more work is needed to show that it can be solved in a full 
neighborhood of p _1 (0). To do this, we use ideas from [7] and find a smooth 
function m so that instead of (|24p . one has {ma, mb} = in a neighborhood 
of a = b = 0. The flows of the Hamilton vector fields H ma and Pl m b will then 
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commute, and one can compute a change of coordinates in a neighborhood 
of p -1 (0) which reduces H mp to d yi + id V2 . After this, the equation (fT8|) is 
easy to solve. The details are given in [16] . 

6. Partial data results 

Let Q, C R n , n > 3, be a bounded simply connected C°° domain with 
connected boundary. We want to discuss the unique determination of the 
magnetic field dA and electric potential V, if one knows the boundary mea- 
surements Aa,v on a subset of the boundary <9S7. 

Let xq G R n \ ch(f2), where ch(fi) is the convex hull of Q. We define the 
front face of d£l relative to xq by 

F(x ) = {xedn-. (x- xq) ■ v{x) < 0}, 

and we take F to be an open neighborhood of F(xo) in Oft. Also, let A ta , n 
be the tangential component defined in 

The following result was proved for C 2 magnetic potentials in [6], and it 
was extended to Holder continuous potentials in [16] . 

Theorem 6.1. Let Aj G C £ (Tt;C n ), e > 0, and Vj G for j = 1,2. 

Also assume that is not a Dirichlet eigenvalue of H A .y. in $7. If 

A Al yJ\ p = A Aa ,v a f\p for all / G H 1 / 2 ^), 

then dA\ = dA2 in O and V\ = V2 in f2. Also, (Ai) tan = (^42)tan ° n 9Q. 

We describe the proof in the case where Aj G C 2 (ft; C n ) with Aj ■ v = 
on <9S7, following [6]. See [16] for the C e case and some further details. We 
start by choosing u\ and 112 to be CGO solutions, given by Theorem 15. 1] to 
the equations H^yui = and H^ 2 y 2 U2 = 0. They will be of the form 

ui = eK(v+^)(oi + n), 
U2 = eK(-^)(a 2 + r 2 ), 

where ip is the logarithmic weight (I22p . ?/> is the distance function ([31H . and 
are amplitudes solving transport equations. One reason for working with 

limiting Carleman weights is that u\ and 112 can be chosen with opposite 

signs in front of ip. With this choice, the terms e^h^ which grow as h — > 

cancel in the product uiU2- 

Take 112 to be the solution to Ha^Vi^ = with U2 = ui on dfl, and let 

u = u\ — U2- Integration by parts shows 

(H A2 y 2 u\u 2 ) = -(d u u\u 2 )dn. (33) 

Since the DN maps coincide on F, and since d u u = (A^v^ — A A2 y 2 )m, one 
sees that d v u = on F. Thus the right hand side of ([33l) reduces to an 
integral over <9f2 \ F, and one would like to estimate d v u in this set. 

As shown in [3], one can estimate d v u by using a Carleman estimate with 
boundary terms. This is an extension of (|26p to functions v which vanish 
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on dft, but are not necessarily compactly supported in £1. The estimate has 
the form 

Vh\\ y/d u ipe~h d u v\\ L 2( dn+ ) + ||e~fcv|| + ||e~fc/iVu[| 

< Ch\\e~hH A2 y 2 v\\ + CVh\\y / -d 1/ tpe~hd 1/ v\\ L 2( dn _y 

Here d£l± = {x £ d£l : db d u ip(x) > 0} are the front and back sides 
of d£l. This estimate can be applied to u, and one gets a bound for 
h d v u\\ L2 , gQ ^py Consequently 

h(d u u\u 2 )dn -> 

as h — > 0. Computing the limit of h times the left hand side of (I33p . and 
using the explicit forms for u\ and u 2 , yields 

\/(cp + iip) ■ (Ai- A 2 )adx = 0, (34) 

Jn 

where a = liuih^Q a\d 2 . The expression on the left may be thought of as a 
counterpart to (I20p . and we will see below that (I34p asserts the vanishing 
of a certain nonlinear Radon transform of d(A\ — A%). This will imply 
dAi = dA 2 . 

We switch to the complex notation used in Section [5l Thus, z = x\ + ir 
where x = (xi,r9). We also write P$ for the two-plane consisting of points 
(xi, rO) with 9 fixed, and fig = Q D Pg. The formula (|34|) becomes 

/ (l -(ei+»er) • (Ai- A 2 )(z,9)e iiS> dz Adz) d9 = 0. (35) 
Js n ~ 2 ^ Jn g z ' 

Here <!> arises from the amplitudes a\ and 02, and will satisfy 

dg$(z,9) = -i(ei + ie r ) • (vli - A 2 )(z,fl). (36) 

We need a slightly more general version of (|35p . From the transport equation 
for 01, one sees that a\ can be replaced by a\g\ where dzgi(z, 9) = 0. Then 
(|35l) holds with e** replaced with e l ^g\. Choosing gi(z,9) = zg{z)g{9) with 
$2<? = 0, and varying g, we see that for almost every 9 one has 

f ( ei +ie r ) ■ (A^ A 2 )(z,9)e i ^ z ^g(z)dz Adz = 0. (37) 
Jn g 

We will show that (|37p remains true with e J * replaced with 1, which 
corresponds to converting a nonlinear Radon transform to the usual one. 
This will be done by " taking a holomorphic logarithm" . The procedure has 
similarities with [29] . The first step is to use the equation (|36p and integrate 
by parts in ([37]) . which yields 

I e i0 g(z) dz = (38) 

Jdflg 

for g holomorphic in Qg. This orthogonality condition means that e'*|gn e 
is the boundary value of some holomorphic function F E C(Qg). It can be 
shown that F is nonvanishing and has a holomorphic logarithm G G C(Q,g). 
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Taking g = % in (|38j) . and noting that G = i$ on <9Q up to constant, one 
gets 

/ $ dz = 0. 

Integrating by parts and using ([36]) again gives 

(ei + ie r ) • (^i - vl 2 )(2, 0) dz Adz = 0. 

Going back to x-coordinates, and varying xq and u slightly, this will imply 

( C-(Ai- A 2 ) dS = (39) 
Jpnn 

for all two-planes P with d((0, ei), T(P)) < <5, and all £ G P. 

Finally, we indicate how (j39[) implies dA\ = dA 2 in the case n = 3. The 
argument uses analytic microlocal analysis. By approximation, one may 
assume Aj £ C7?°(R 3 ) in ([39]). Fixing £,rj £ R 3 and writing A = A\- ^4 2 , 
we obtain from ([39]) 

y (dA^Arj) dS = 

for all P with d((0, ei), T(P)) < 5. By the microlocal Helgason support 
theorem, the vanishing of the Radon transform at the planes P implies the 
microlocal analyticity of (dA, £ A rj) on the conormal bundle of every such 
P. Consequently, if {dA^Arj) is nonzero, then it is microlocally analytic at 
some point of the normal set of supp((cL4,£ A r/)). But this contradicts the 
microlocal Holmgren theorem, so one must have dA\ = dA 2 . We refer to [6] 
for more details, and for the case of electric potentials. 
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